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Abstract 


The amplitude equation for an unstable electrostatic wave is analyzed using 
an expansion in the mode amplitude A{t). In the limit of weak instability, i.e. 
7 ^ 0 '*' where 7 is the linear growth rate, the nonlinear coefficients are sin¬ 
gular and their singularities predict the dependence of A{t) on 7 . Generically 
the scaling |A(t)| = 7 ^/^r( 7 t) as 7 —> 0 ^ is required to cancel the coeffi¬ 
cient singularities to all orders. This result predicts the electric field scaling 
|Afc| ~ 7 ^/^ will hold universally for these instabilities (including beam-plasma 
and two-stream configurations) throughout the dynamical evolution and in 
the time-asymptotic state. In exceptional cases, such as infinitely massive 
ions, the coefficients are less singular and the more familiar trapping scaling 
\Ek\ ~ 7 ^ is recovered. 
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I. INTRODUCTION 


A. Overview and background 

We recently presented a detailed analysis of the amplitude equation for an unstable 
electrostatic mode in an unmagnetized Vlasov plasma (henceforth (I)). [Q] In that work, the 
amplitude equation was studied using an expansion in the wave amplitude, and only the 
electron motion was considered with the ions treated as a hxed neutralizing background. In 
this paper we study the effect of including mobile ions and investigate the mode amplitude 
equation for an electrostatic wave in a multi-species unmagnetized plasma. As in (I), we 
focus on the regime of weak instability and carefully analyze the structure of the amplitude 
equation in the the limit 7 ^ O’*" where 7 is the linear growth rate for the unstable wave. 

Incorporating the effects of hnite ion inertia is a signihcant step for several reasons. With 
hnite mass, the ion motion occurs on a hnite timescale Tj < cxo and in the limit of very weak 
growth rates this will be a fast timescale, i.e. T* << 7“^ as 7 ^ O’*". Over the duration of the 
instability, the ions can readily respond to the electric held of the wave especially resonant 
ions near the phase velocity. This ion response is always present even for high frequency 
modes such as plasma waves since at the phase velocity the wave frequency is Doppler-shifted 
to zero. In addition, there are modes such as ion-acoustic waves that depend on ion density 
oscillations and are suppressed when the ions are treated as hxed. Allowing for mobile ions 
is a prerequisite, if we seek the amplitude equation for an unstable ion-acoustic mode. 

This paper can be read independently, but the reader is referred to (I) for a discussion 
of the history of the problem and additional introductory material. An essential diherence 
between our approach and previous work lies in the choice of unperturbed state. Earlier 
theories assumed an equilibrium with a neutrally stable mode and obtained ill-dehned ex¬ 
pansion coefficients. p|- 0 This can be avoided by taking the weakly unstable equilibrium as 
the unperturbed state; a choice that naturally leads one to work with the unstable manifold. 
A brief summary of our results on the multi-species case has appeared elsewhere. 0 

In the limit 7 —>■ O’*", the amplitude equation dehnes a kind of singular perturbation 
problem whose detailed features reveal asymptotic scaling behavior of the nonlinear wave. 
This is a key idea behind our approach and it can be formulated more precisely as follows. 
The mode eigenvalue A = 7 — za; can be complex (beam-plasma) or real (two-stream), but 
in either case the equations for the amplitude A{t) = p{t) have the form 

p = 'jp + aip^ + a2p^ + O(p^) (1) 

e = uj + a[p^ + a'^p^ + 0{p^). (2) 

Since the Vlasov equilibrium is assumed to be spatially homogeneous, the evolution of p{t) 
decouples from the phase 9(t) . Suppose for simplicity that the coefficients aj are constant 
(ignoring their dependence on 7), then we may view (J^) as a singular perturbation problem 
with the linear term 7p representing the perturbation. If 7 > 0 , then there is always a 
neighborhood of the equilibrium p = 0 where the perturbation dominates the unperturbed 
system and can completely change the dynamics. 

As in other such singular problems, a possible strategy is to seek a (singular) change of 
variables which transforms m into a regular perturbation problem. Thus it is natural to 
dehne a new amplitude r(r) by 
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pit) = 7^(7'5t) 


(3) 


and rewrite (|^ in terms of r(r) 

^ = 7i-'5r(r) + 72 / 3-501 (7)r3 + 7^/3-%(7)r5 + • • ■. (4) 

dr 

If possible, the choice of [3 and 5 shonld be made so that (z) each term in (|) is well 
behaved as 7 —O'*" and (m) the effect of 7 > 0 is a regnlar pertnrbation of the system at 
7 = 0. The latter condition reqnires 5 = 1; otherwise the linear term will continne to be 
a singnlar pertnrbation. The exponent (3 then needs to be chosen (if possible) to achieve a 
balance between the nonlinear terms and the linear term. In the simplest sitnation, when 
the coefficients Oi, 02 ,... have well-dehned finite limits as 7 — 0 , then one has the exponent 
(3 = 1/2 and this yields 

(IT 

— = r(r) + ai( 7 )r^ + 7 a 2 ( 7 )r® H-. (5) 

CLT 

The nnpertnrbed system f = r + ai(0)r^ now inclndes the linear term and near r = 0 the 
effect of small 7 is nnimportant. The change of variables (|]) is singnlar at 7 = 0 as expected, 
bnt nevertheless we have determined the overall scaling (3 = 1/2 for nonlinear solntions near 
the eqnilibrinm in the regime of weak instability. In addition, since terms in (|^) at hfth 
order or higher are at least of order 0 ( 7 ), the nnpertnrbed problem trnncates to a simple 
balance between the linear instability and the dominant nonlinear term. 

For an nnstable electrostatic wave, the sitnation is qnite different becanse the nonlinear 
coefficients blow np as 7 —>• 0 and the exponent (3 = 1/2 does not yield a well-behaved 
eqnation for r(r). We attack the problem of hnding a satisfactory valne for (3 in two steps. 
First, we determine what choice for (3 will control the singularity of the cubic coefficient ai 
in (|l|), and then we investigate whether this choice will also remove the singularities of the 
higher order nonlinear terms. Typically, the cubic coefficient diverges asymptotically like 
7 “^ in which case we must choose (3 > 5/2 in order to obtain a hnite cubic term in the 
rescaled amplitude equation (^. 

Our analysis shows that setting (3 = 5/2 yields a theory that is hnite to all orders 
as 7 —0 at 7 = 0 retains (rescaled) nonlinear terms that balance the linear instability. 
Choosing a larger exponent (3 = b/2 + e would also yield a hnite theory to all orders, but 
now the cubic term is of order 0 ( 7 ^'^) and all the higher order terms (j > 1 ) are of order 
there is no balance between the linear and nonlinear terms at 7 = 0. For e > 0, the 
larger exponent suppresses nonlinear ehects too strongly. In this sense, the scaling 

pit) = 7/5/V(7t) ( 6 ) 

is uniquely determined. Of the two steps, the second is by far the most difficult; it is relatively 
straightforward to compute ai and determine its divergence (Section |In[). Verifying that a 
given choice of (3 works to all orders requires a detailed study of the recursion relations of 
the theory and the hnal theorem is proved by induction (Section 0 )- 

In special cases, the cubic coefficient is less singular, diverging like 7 “^ instead of 7 “^. 
The limit of inhnitely massive ions is the best known example of this situation, and in (I) we 
found that (3 = 2 was the scaling needed to make the nonlinear amplitude equation hnite to 
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all orders. In this paper, we identify two additional examples where the cubic coefficient has 
a 7 “^ singularity: if the effects of the resonant ions are negligible (in a sense we make precise 
in Section |111 Dj) or if the electron and ion masses are equal, i.e. for an electron-positron 
plasma. In the former case we prove that the exponent /3 = 2 is correct, but the precise 
value of (3 for the electron-positron system is not unambiguously determined due special 
complications in the singularity structure of the amplitude expansion. 


B. Notation 


Our model is a one-dimensional, multi-species Vlasov plasma dehned by 


dt 


+ V 


dx 


+ E 


dv 


0 


(7) 


dE 

dx 


f dv E^^\x,v,t). 

a J — OO 


( 8 ) 


Here x, t and v are measured in units of u/uje, and u, respectively, where m is a 
chosen velocity scale and cUg = Ane^ne/me- The plasma length is L with periodic boundary 
conditions and we adopt the normalization 



(9) 


where Qs = eZg is the charge of species s and = gsme/emg. Note that = —1 for 
electrons and that the normalization (^) for negative species makes the distribution function 
negative. 

Let Eo{v) and f{x,v,t) denote the multi-component helds for the equilibrium and per¬ 
turbation respectively and k the matrix of mass ratios. 



/ /(«!) \ 
f{s,) 

Eo = 

Fli“' ) 

p(.,) 

K = 


[ ■■ j 


\ J 

V 


0 0 ■■■\ 
0 0 ... 


V : 




then the system (|^) - (|p can be concisely expressed as 

| = £/ + V(/) 


where the linear operator is dehned by 

OO 

C/= E e“nL,f,){v 

/ = —OO 

(Limv) = ^ 


-il 


0 / = 0 
ds'), 


Vfl{v) + K ■ T]i{v) Es' f-oo dv' fl \v') I ^ 0, 


( 10 ) 


( 11 ) 


( 12 ) 

(13) 
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with rii{v) = —dyFo/l"^, and the nonlinear operator Af is 


mf)= T. 


T. 


l=—oo 


K 


dfm-l 

dv 


Y. r dv'ft'\v'). 

/ J—oo 


( 14 ) 


In the spatial Fonrier expansion (|^, I denotes an integer mnltiple of the basic wavennmber 
27r/L, and a primed snmmation as in (p!4D omits the I = 0 term. The notation k, ■ rjiiv) or 
K ■ dyfm-i denotes matrix mnltiplication. 

An eqnilibrinm Fo{v) depends on parameters snch as the densities or temperatnres of 
the varions species, bnt this dependence will be snppressed nnless it is explicitly needed. 

Symmetries of the model (|^) and the eqnilibrinm Fq{v) are important qnalitative 
featnres of the problem, and play a significant role in Section || when we formnlate 
the amplitnde expansions. Spatial translation, Ta : {x,v) {x + a, v), and reflection, 
TZ : {x,v) (—X, —v), act as operators on f{x,v,t) in the nsnal way: if a denotes an arbi¬ 
trary transformation then (a f){x,v) = {.x,v)) where (a f){x,v) denotes the trans¬ 

formed distribntion fnnction. The operators C and Af commnte with 7^ dne to the spatial 
homogeneity of Fq, and if Fq{v) = Fq{—v), then C and Af also commnte with the reflection 
operator IZ. With periodic bonndary conditions, x is a periodic coordinate so Ta and IZ gen¬ 
erate 0(2), the symmetry gronp of the circle. If only the translation symmetry is present, 
then the gronp is SO(2). 

An inner prodnct is needed in Section || to derive the amplitnde eqnation. For two mnlti- 
component helds of (x, n), e.g. B = ..) and D = 

we dehne their inner prodnct by 


/ L/2 1*00 

dx / 

^ -L/2 J—oo 


rLl2 

dvB^^\x,vyD^^\x,v) = dx <B,D> 

J-L/2 


(16) 


where 

/ OO 

dvB^^\x,vyD^^\x,v). 

-OO 


(16) 


C. Summary of linear theory 


The spectral theory for C is well established, and the needed resnlts are simply recalled 

Tod The eigenvalnes A = —ilz of C are determined by the roots 


to establish onr notation. 

Ai{z) = 0 of the “spectral fnnction 


A,(.).i+f,.4AAhM, 

J-oo V — Z 


(17) 


If the contonr in ([T7|) is replaced by the Landan contonr for Im(^) < 0 then we have the linear 
dielectric ei{z); for Im(^) > 0, Ai{z) and ei{z) are the same fnnction. For SO(2) symmetry, 
the eigenvalnes are generically complex; for 0(2) symmetry, when the eqnilibrinm is also 
reflection-symmetric, then real eigenvalnes can occnr as in a two-stream instability; see for 
example fH . 
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Associated with an eigenvalue A = —ilz is the multi-component eigenfunction = 




^ • m 

V — z 


( 18 ) 


There is also an associated adjoint eigenfunction T(a;, n) = e^’'^'ilj{v)/L satisfying (T, T) = 1 
with 




1 

A'i{z)*{v — z*) 


(19) 


Note that all components of ‘ip{v) are the same. The normalization in (|T^ assumes that 
the root of Ai{z) is simple and is chosen so that Kip,ip >= 1. The adjoint determines the 
projection of f{x,v,t) onto the eigenvector, and this projection dehnes the time-dependent 
amplitude of T, i.e. A{t) = (T, /). 


II. AMPLITUDE EQUATION ON THE UNSTABLE MANIFOLD 

A. Unstable linear modes 

The equilibrium Fq{v) is assumed to support a “single” unstable mode in the follow¬ 
ing sense. We shall assume that E^, the unstable subspace for C, is two-dimensional. With 
translation symmetry and periodic boundary conditions, this is the simplest instability prob¬ 
lem that can be posed. 

Henceforth, let k denote the wavenumber of this unstable mode that is associated with 
the root Afc(^o) = 0 which we assume to be simple, i.e. A'j^{zq) ^ 0. The corresponding 
eigenvector is 

= ( 20 ) 

\ V — ZqJ 

The root zq = Vp + i'y/k determines the phase velocity Vp = ui/k and the growth rate 7 of 
the linear mode as the real and imaginary parts of the eigenvalue A = —ikzo = 7 —io;. From 
(Q) and Ak{zo) = 0 , we obtain 

Aiizo) = {l^ - A)^)//^ (21) 

this identity will be useful shortly. 

When Fq{v) lacks reflection symmetry, then Ti typically has a non-zero phase velocity 
and A is complex. In this case, the identities Ak{z) = A_k{z) and Ak{z)* = Ak{z*) imply 
three additional modes: T 2 , and T 2 where 

<i> 2 {x,v) = F^^'iP{vy. (22) 

These eigenfunctions correspond to eigenvalues A*, —A*, and —A, respectively, and £11 out 
the eigenvalue quartet characteristic of Hamiltonian systems. In the absence of reflection 
symmetry, the eigenvalues are typically simple and the unstable subspace is two-dimensional. 
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When Fq{v) is reflection-symmetric, the eigenvalues may be either real or complex. In 
either case tl'j and TZ ■ d'j are linearly independent so the eigenvalues have multiplicity two. 
Now is two-dimensional only in the case of a real eigenvalue since a multiplicity-two 
complex conjugate pair implies a four-dimensional unstable subspace; this latter possibility 
is not considered further. 

The components of the distribution function along the unstable eigenvectors Ti and 
are identified by writing 


/(x, V, t) = [yl(t)\['(x, v) -f cc] -f S{x, V, t) 


(23) 


where A{t) = (T, /) is the mode amplitude for T and (T, S) = 0. In (pdf ), the subscript on 
Ti has been dropped, and T = exp{ikx) 'll)/L is the adjoint function for zq from (H)- 

The action of translation 7^ and reflection 7Z on f{x,v,t) implies an action by these 
operators on the variables (24,5'). From (^3|) we note that Taf{x,v,t) = f{x — a,v,t) is 
equivalent to 


T, (A, S{x, v)) = S{x - a, v)), (24) 

and reflection 7lf{x,v,t) = —v,t) yields 

TZ {A, S{x,v)) = {A*, S{—x,—v)). (25) 

When these transformations are symmetries, these relations are important for organizing 
the amplitude expansions introduced in Section |n . 


B. Derivation of the amplitude equation 


In the (A, S) variables, the Vlasov equation (^) becomes: 


A 

dS 

'm 


XA + i^,Afif)) 

CS + Afif)-\i^,Afif))^ + cc 


where 





__ roo , 

/ J-oc 


(26) 

(27) 


(28) 


In writing (|^) we have used the adjoint relationship (T,£S') = (£lT, ^) = A(T,S') = 0 
and in (1^) an integration by parts shifts the velocity derivative onto ^fJ. These coupled 
equations are equivalent to (pdj) ; however by restricting them to the unstable manifold we 
obtain an autonomous equation for 24(f). This procedure is briefly summarized; a more 
detailed description is provided in (I). 

The unstable subspace is invariant under the linear dynamics dtf = Cf, but this 
no longer holds once the nonlinear terms M{f) are included. Instead we assume there is a 
two-dimensional invariant surface, the unstable manifold that represents the nonlinear 
deformation of E^. The unstable manifold is tangent to E^ at the equilibrium, and near 
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Ff) it can be described by a function F[{x,v, A, A*) which locates the manifold “above” the 
point in the geometry is illustrated in Fig. 1. 



Fig 1: Schematic representation of the unstable manifold 

Thus 


P{x,v,t) = [A{t)'^{x,v) + cc] + E[{x,v,A{t),A*{t)) 


(29) 


represents a distribution function on PF’' 
mined from H, i.e. 


In this expression, the evolution of S be deter- 


S{x, V, t) = H{x, V, A{t), A*{t)) = 


f H^^^^{x,v,A{t),A*{t)) \ 
H^^^\x,v,A\t),A*\t)) 


V 


(30) 


/ 


When this representation is substituted into (^) - (pTl) we obtain 


~dt 


A = AA + (^,Ar(r)) 

= £H + Af(r) - + CC 


(31) 

(32) 


Note that (|^) dehnes an autonomous two-dimensional flow describing the self-consistent 
nonlinear evolution of the unstable mode; this is the amplitude equation we shall study. 
Certain general features of the amplitude equation follow from the transformations (^41) 
when these are symmetries of the Vlasov equation (PTj). For our problem, the am¬ 


plitude equation always has translation symmetry (^^ and we can apply standard re¬ 
sults on the form of such symmetric equations. In particular, a two-dimensional vector 
held A = V{A,A*) that is symmetric with respect to ^4 —>■ can be written as 

V{A,A*) = Ap{a) where a = \A\^ and p(cr) is a smooth function determined from V. K 
Hence we know the right hand side of ( 0 ) takes the general form. 


Ap(a) = \A+(9,Ar(r)), 


(33) 


where the function p(a) must be determined from the Vlasov equation. Typically p(a) is 
complex-valued, however when Fq is rehection-symmetric then p(a) is forced to be real. 
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C. Analysis of H{x,v, A, A*) 


An equation for H follows by requiring consistency between 
time derivative of 


]) and 


with the right hand side of (|3^ gives 


dH 


dH 


A+^^A* =CH + M{n - [(vb, M{n) ^ + cc 


Equating the 


(34) 


which is to be solved for H subject to (^). Since the manifold is tangent to at the 
equilibrium, H must satisfy 


0 = n, 0, 0) 


OH 


x, n, 0, 0) 


OH 

dA* 


X, V, 0, 0). 


(35) 


The symmetries of the problem impose general restrictions on H because the unstable 
manifold is mapped onto itself by a symmetry transformation. For a translation this in¬ 
variance means that /“(x — a,v,t) corresponds to a point on the manifold and thus can be 
written as, 


(Zi/“) {x, v) = [AT(a; — a, v) + cc] -|- H{x — a, v, A, A*) (36) 

= [A'T(a;, v) + cc] + H{x, v, A', A'*). (37) 

Consistency between these two expressions requires that the transformed amplitude A' is 
related to the initial amplitude Ahy A' = e~^^°‘A (cf. (|^) ), and this implies that H satishes 

H{x - a, n. A, A*) = H(x, v, 6"*^“A*). (38) 

If, in addition, the system (^) - (^) has reflection symmetry, then similar reasoning yields 

H{-x,-v,A,A*) = H{x,v,A\A). (39) 


These general symmetry properties of H usefully constrain the form of the Fourier ex¬ 
pansion 


OO 

H(x,v,A,A*)= e^'^Hi(v,A,A*) 

l=—oo 


l=—oo 


H^"\v,A,A-) 


m 


Applying (|3B|) shows that the components Hi must vanish unless I is an integer multiple of 
k and that the non-zero components have the general form 


Ho(v,A,A*) = aho(v,a) 

Hk(v, A, A*) = Act hi(n, a) (41) 

H^k{v, A, A*) = A”" hm{v, a) for m>2 

where H^i = Hi*. The functions hm are not determined by symmetry; however, if reflection 
symmetry also holds, then 


hm(-v,a) = hm(v,a)*. 


(42) 
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The reasoning leading to these symmetry results can be briefly outlined. The translation 
symmetry (^) requires 


A, A*) = Hi{v, (43) 

By choosing a = 27r//c, this implies we And Hi = 0 unless exp(27ri(//fc)) = 1, hence the 
non-zero components satisfy I = mk for integer m. Now recall that a function G(n, A, A*) 
which is invariant under A —>■ Ae~^^°' has the form G{v, A, A*) = g{v, a) where a = \A\^ and 
g{v, a) is a smooth function determined by G. The relation in (^3D implies that G{v, A, A*) = 
{A*)‘^Hmk{v, A, A*) is such a invariant function. The factor (A*)'" in G requires that the 
corresponding g{v,a) contain an overall factor of cr™; hence we set g{v,a) = a'^h{y,a). 
The resulting expression, (A*)^Hmk = <j"^h{v,a), then implies Hmk = A'^h{v,a). This 
general form, together with the requirement (|55D, leads to general expressions for the Fourier 
components of H in (p|). 


D. Equations for p{a) and hm{v,a) 


The calculation of the amplitude equation (|31|) reduces to a determination of the function 
p{a) in (p3D; this requires the evaluation of (T, A/'(/“)) from (pS]) and involves the Fourier 
components of /“ 


/“ = [Ai){v) 6i^k + A*^/J{v)* 6i_k] + Hi{v, A, A*). 
Combining (p8|) with (|^) and (^4|) yields 


(T,Ar(r)) = 


-iAa 


F 2 

< K- (ho - ^ 2 ) > +— < K-'ip* > 


(44) 


(45) 


+a 


Fi < dyip, K- ho> —F)" < dyip, K ■ h2> +— < dyip, K ■ h\> 


< Byp), K ■ /l3 > 


00 ^1-2 


+ ^ ■ ^*1-1 > ^ ■ ^'+1 >] 


1=3 


In this expression, sum over the velocity integrals of is denoted by 

/ OO 

dv h^^(v,a). 

-00 

Now comparing (|3^) and (^5]) provides the desired expression for p 

%CT f r 

p((j) = A —r 1 < (^0 - ^ 2 ) > < Bytp, K-ip* > 


(46) 


(47) 


+a 


k 

r ~ r* ~ 

Fi < dyip, K- ho> —F); < dyip, K- h2> +— < dyip^ K- h\> —^ < dyip^ K ■ ho > 


00 ^1-2 

1=3 ^ 
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This expression for p involves the functions hm{v,(j) which define H. The equations 
that determine hm{v, a) follow from the general equation for H in (^4|) . In presenting these 
relations we shall use the orthogonal projection defined by the unstable mode. Let V 
denote the projection operator onto 'ip{v) of a multi-component function g{v), 

{Vg){v) =<'^,g> ijjiv), (48) 


with the orthogonal projection V± = I — V given by 

CP± 9 )(«) = g(.v)-ip(v) <i>,g> . 


(49) 


With the previous expression for (T,A/'(/“)) in (^) and the notation in (^) for the 
Fourier components of if, the components of take the form: 


{p + p* 


ho (j 


dho 

da 


i d 
k dv 


K ■ 


oo 1-1 

r + - ipri) + f/Ti + ^ 

1=2 ^ 


— cc, 


(60) 


(2p + p’)hi - (Lthi) + (p + p")ct 


dhi 


da 


^ ■ i ^0 ~ ^2 + 2 '^* r 2 + (T 


(51) 


fobi - hTi + -h*T 2 - -hr; 


00 0-1—2 

E -r - ‘"'“+1^' 

1=3 ^ 


2p h 2 - L 2 kh 2 + {p + p*)cr 


dh2 


i d 
k dv 


K ■ lijj + a 


hi + tliTi- h3 + ^hoT2 


a 


fiTi - fgrj - 


„ 00 1-2 

-=hn + E -r 

l=i ^ 


(52) 


and 


mp hm - Lmkhm + (f» + P*)(T 


dh^ 


id J ^ , V’ T. , hm-i^ 

-K ■ \ hm-1 H-7! m-1 + 2^ —^ f I 


k dv 


m — 1 


1=2 


(53) 
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+a 




h r h , r , 1 r 

1 ~ n^m+l H-7-1- m-1 H-J- m H-—T-L m+1 

m — 1 m m + 1 


-hm+l^\ + 


w 


m + 1 


m+1 


+ E 


oo ^l—m 


htrJl 


l=7n-\-2 


oo I 

- E 

1=2 ^ 


for m = 0,1, 2, and m > 2, respectively. 

Together with 0) these component eqnations determine the functions p(cr) and 
{hm (o') 15 ^= 0 ) however they cannot be solved except using the expansions introduced in 
the next section. From a practical standpoint, we have achieved a reduction of the problem 
to the analysis of functions of a single real variable, i.e. a. In the study of the amplitude 
equation (|^) this reduction represents a very useful simplihcation. 


III. EXPANSIONS, RECURSION RELATIONS, AND PINCHING 

SINGULARITIES 


The amplitude equation on the unstable manifold, 

A = Ap{a), (54) 

is analyzed by expressing p(o'), hm(n, a), and rm(o') as power series in a: 

CCi OO OO 

P(^) = E PJ ^) = E r„i(cr) = ^ Tm,j cr^ (55) 

j=0 j=0 j=0 

where = f dvEs ^mj (^)- 

The coefficients pj and h^giv) are calculated by recursively solving (^71) and (^0]) - (p^. 
A recursion relation for pj follows from (^^ by inserting the expansions above and solving 
at each order in a; this calculation yields 


Po = X 

for i > 1 


( 66 ) 

(57) 


where 


Aj 


B,= 


■ < K ■ (hoj-1 - haj-i) > “2 ^ K-^J* > Taj-i 

t-2 

“ E [< ■ ho,j-i-2 > Ti,/- < dy'ip, n ■ h2,j-i-2 > Tl, 

1=0 

E h < Svi’, K ■ h’id-i-2 > r2,i 


(58) 


1=0 


(59) 


E 

m=0 


F F* 

^ j-l+l,m ,07 1 * , j-l,m 


J - l + l 


< a+r « ■ h*_i > - > 

J - L 
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Here and below, a summation is understood to be omitted if the lower limit exceeds the 
upper limit. This organization of the terms in (^) and (^) will turn out to distinguish to 
different singular behaviors: Aj ~ ^ ^-(5i-2)_ 

The corresponding relations for hm,j{v) are determined similarly by applying the expan- 


-1 


sions (p5D to the general equations (|^ - (|^); the resulting equations are solved for hm,j{v) 
Our results are conveniently stated in terms of the resolvent operator Ri{w) = {w — Li) 
and certain auxiliary functions Im,j{v) whose detailed expressions are provided below. For 
m = 0, hoj is simply 


ho,j{v) = 




(l+j)(A + A*)’ 


and the coefficients hm,j for m > 0 require the resolvent 


RmkiWfnj) I’i 




(60) 


(61) 


with the complex numbers Wm,j defined by 

Wm,j = (j + ^m,l)(A + A*) +m\ = 2{j + 6m,l) 7 + "iA. 


(62) 


A general expression for Riiw) follows from (JT^ by solving (w — R) f = g for f. 

Here both g{v) = ,...) and / are multi-component helds, and Ri{w) acts by 


Ri{w)g = 


( {Ri{w) g)^^^\v) \ 
{Ri{w) g)^^^\v) 


V 


(63) 


/ 


where 


{Ri{w) g)^"\v) = - 


il{v — iw/l) 


gA(v)- ^ y 

Riilw/l) ^. 


dv' 


v' — iw/l 


(64) 


Note the following signihcant feature: the arguments Wm,j of the resolvent in ([6l|) deter¬ 
mine poles of hm,j{v) located at v = Zm,j where Zm,j = iwm,j/fnk. For m > 1, these poles 
always fall in the upper half-plane above the phase velocity Vp: 




— Zq + 


Rd. 


m,j 


= V„ + 


*7(1 -|- d. 


m,j, 


k 


(65) 


where dmj = 2(j + S„a)/m. 


A. Recursion relations for L 


m,] 


The auxiliary functions Io,j{v) are dehned by 

hoiy) = 


( 66 ) 
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j = 0, and for j > 1 


hA'^) = - II (1 + n)iPj-n + P*j-n) ho,nA) 

n=0 

i-2 


i d 


Kj-l - '/’r'lj-l + ^ K,n^i,j-n-2 


n=0 


j-1 n U* 

h>h 3-n +I 


r,_ 


ji—n+l,n—/ 


— cc 


( 67 ) 


For m = 1, the functions are given by 




= - H [(2 + + (1 + n)p*_n] hi 


n=0 




2,j 


r 1 

+ H ho^nhlj-n-l — h2^rJ'ij_n-l + '^hi,^T2,j-n-l 

n=0 ^ 

1 ^ / h* \ 

i-2 n 

-EE 


rF. 


j—n-\-2^n—m 


n=0 m=0 _ 


hj-n-\-2,m 

j - n + 1 


( 68 ) 


The notation V^giy) was dehned in (|49|): V^giy) = g{y) — 'ip{v) Kip^g >. 
For m = 2 the functions l 2 ,j{v) are given by 




(69) 


for j = 0, and for j > 1 


i-i 


h,j{v) = - E [(2 + APj-n + np 


U h 


n=0 


2,n 


i d 


i-i 


+ ^ ■ ■{ ^l,i-l + — hsj-i + —'0*F3 j_i + — ^ ho,nr2j-n-l 


n=0 


(70) 


t “2 r 

+ E/ ^linhl j_n,-2 ~ — —/i 4 „F2 j_„_2 + 77^4 ,^F3 j_„_2 

n=0 L / O 


i-3 

E- 

n=0 


^5,n 


m=0 

i -4 n 7 

■p* \ \ 'v —77 + 2,772 -p:): 

3,j-n-3 — ^ 


h* 

j—n,m -p 

, . I 7—nH-2,n—m 

J - n + 2 


E- 


n=0 777=0 ^ ^ 


j—n^n—m 
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( 71 ) 


For m > 2, the functions Im,j{v) are given by 




+ n)pj_n + np*_^ hm,n{y) 


n=0 


m-2 j L 

' —/. 


ip* 


l,j-n 


hm+l,j — l “1“ ^ hm-i-ij — i 


i-1 

E 

n=0 . 
i-2 

E 

n=0 


m + 1 

7 -p I ^l,n p . ^0,71-p 

l,n-*- l,j—n—l “i“ m—l,j—n—l “r i 77^ 72—1 

m — 1 m 


h* 

U "P* I P 

"^m+l,nJ- 1j_7x_2 “t | m+lj—n—2 


m + 1 


n / h* h 

'^j—n,l p '^m+j—n,/p* 

-L j-\-m—n,n—l 1 -L 


E I ■ 

\j + m-n 


j - n 


j—n,n—l 


in this last expression, if a subscript is negative the term is understood to be omitted, e.g. 
for j = 0, hm+i,j-i is omitted. 


B. Useful identities 


The relations (0. (0) - (0). and (^) - (|7I|) can be applied systematically to calculate 
Pj and hm,j to any order. The leading coefficient po is determined by linear theory (|56|) and 
from the linear eigenfunction ^p one can also calculate ho,o and ^ 2 , 0 , c.f. (pB|) and (§^), re¬ 


spectively. These two coefficients then suffice to calculate pi from dSTj) . From {pi, hop, ^ 2 , 0 }, 
the coefficients hip and hap can be determined and then ho,i and h 2 ,i. This provides the 
input to calculate p 2 , and from this point on the structure of the calculation to all orders 
falls into the simple pattern summarized in Table I. 

For the remainder of the paper, our analysis of the expansion is facilitated by several 


identities which we summarize here. These relations allow F^y, hmj, and certain important 
integrals to be obtained very simply from the auxiliary functions Imj- 

If m = 0 we note that (|60D and (|66D - @) imply Fqj = 0. For m > 0, by substituting 
(1^) into Tm,j = f dv J2s rearranging, we obtain 


^ -i/mk ^ 


foo 


V — z. 


{m > 0). 


(72) 


m,3 


With (ff^ for Fmj, the general form for hm,j in (^) can be re-expressed as 



dm,j{v) 


V — z 


m,j 


( ^ ■ Pmk 
- Zm,j 


F 




(m > 0). 


(73) 


In this way F„jj and are written explicitly in terms of Im,j- 

The coefficients pj depend on integrals of the form < k ■ >, and it is helpful to 

express these integrals directly in terms of Imj also. From (^) we obtain the identity 
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< dyl/j, K ■ hoj > = 


dv 


E.-t'-'/S 


OJ . 


27(1 + j)A'k{zo) 2-00 (v - ZqY ’ 


(74) 


similarly, by integrating ([73|) we find 


7 dyijj^ K ' hrfY^ j 7 


—l 


dv 


T 




mkA'f.{zo) j 7-00 (v - zo)^(v - Zmj) 


m,] 






[v - Zo)^(v - Zm,j) 


(75) 


C. Analysis of the cubic coefficient 


It is instructive at this point to evaluate and examine the cubic coefficient pi in detail. 
This coefficient illustrates the occurrence of pinching singularities and the procedures needed 
to evaluate such singular integrals. 


From we have 


Pi = - 


k [ 


r 


2 0 

< dy^j, K ■ (ho,o - h2fl) > < dyi), K-i)* > 


(76) 


so h^fliy) , h 2 fi{v) and r 2 ,o are needed to calculate pi. From (|60D and (^) we have 

1 d 


hofl{v) = - 


k^ dv 


■ Vk 


{v - Zo){v - z^) 


(77) 


where denotes the square of the mass matrix. From 


^) and (|7^ we obtain 


r2,o — 


and inserting (l6^) into ([73D yields 


2 r 


3/c^ 7-00 {v — ^o)^ 


, /X I fK - dyi>\ F2 ,o f K-Vk 

= ^773 I 7— I - — (— 


(78) 


2k‘^ \ v — zo 


(79) 


Here we have made the substitution p 2 k = hfc/4- 

We are concerned with the form of the various integrals in pi in the asymptotic regime 
of small growth rate 7 —O’*". The behavior of these integrals is determined by any complex 
singularities of the integrand that approach the real velocity axis as 7 ^ O’*". We shall 
always assume that whatever complex analytic structure characterizes Fq there are never 
any singularities that move to the real axis in this limit and therefore the singularities in pk 
can be regarded as irrelevant in analyzing a given integral. 

The singularities of interest are associated with the explicit poles in the denominator of 
a given integrand. For example, F 2 ,o has a third-order pole at v = Zq which approaches the 
real axis from above as 7 —> O’*". In this case, there are no corresponding poles in the lower 
half-plane, and consequently no pinching singularity develops; therefore F 2 ,o has a finite limit 
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as the growth rate goes to zero. Similarly, we see by inspection that h2,o has a third-order 
pole at V = Zq and conseqnently the integrand of < dy'ip, k ■ ^2,0 > has a hfth-order pole at 
V = Zq. However, again there is no pinching singnlarity and this integral also has a hnite 
limit as 7 —>■ O’*". 

The remaining integrals, < dyijj, > and < dy'ip, k ■ ho^ >, are more interesting. The 
simplest example is < dy'ip, n ■ 'ip* > whose integrand has a donble pole at v = zo and a 
simple pole at v = zp,: 


< dy'ip, K ■ 'IP* > = 


, hoo (ll- Zo)^(t> -z;)' 


(80) 


As 7 —> O’*", these poles approach the real axis from above and below; this creates a pinching 
singnlarity and the integral has a singular limit. We analyze this integral in two steps 
that can be applied to all such singular integrals: a partial fraction expansion extracts the 
singularity and the dispersion relation is applied to eliminate integrals over the electron 
distribution function. The motivation for the second step is to incorporate the one relation 
between the different distribution functions and to expose the effect of treating the ions as 
mobile. 

The partial fraction expansion of the integrand yields 


{v - ZoYiv - Zp) 
so that our integral becomes 

^ J-^ {v - ZoYiv - zp) V 27 j 

\ r°o dy . rik)^d 

27 / .-/-oo 



^ dv ^ 


dv 


V- Zn 


V- Zo 


(82) 


tl - Zj tl - Zo 


ik 


2'! j (v - Zo)' 


(81) 


On the right hand side of this expression all the integrals are free of pinching singularities 
and therefore have hnite limits as 7 —O’*'. The singularities of the original integral are now 
explicitly shown in the factors 7 “^ and 7 “^. 

The dispersion relation for the unstable mode Ak(zo) = 0 implies the identities 




V-Zo 
(i> - zo)^ 



(83) 

(84) 


where we have used = _x ^nd the primed sum indicates that the electron term is 
omitted. Substitution of (| 8 ^ - (|^) for the electron integrals in 


yields 



dv ■ 'rjkY^'^ 

{v - zo^iv - zp) 



_ . nOO 

5 J —00 


dvVk^ 

V- Zo 





dv dl 


7) 


-00 (n - zof 


( 85 ) 
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Combining this with (^) we conclude that < dyijj, n-'ip* > typically diverges like 7 “^ in the 
small growth rate limit. There are exceptions however; for example, in the limit of hxed 
ions, —>■ 0 for s 7 ^ e, then we are left with 


E 


dv 


'-oc {v - zoYiv - z^) \27 


AU^o). 


( 86 ) 


and the divergence is merely 7 “^. 

We treat the remaining integral < d^'ip, k ■ h^ Q > in the same fashion. This integrand 
now has a fourth order pole at u = 2:0 and a simple pole a.t v = Zq\ 


< K ■ hofi >= 


■E 


dv 


s (n - ZoYiy - z^)' 

Again we make a partial fraction expansion 

1 / ifc E 1 1 

(v-ZoY(v-z^) [2jJ _v-z^ v-Zo_ 

1 (ik 


(87) 


( 88 ) 


ik 


1 


ik 


\27/ {v-ZqY \2'y J {v-ZoY \‘^1 J (v - ZqY' 

and eliminate the electron integrals to hnd 


E 


dv {k^ ■ r]kY''^^ 


ikV 


+ 2 i^'E)(l - E)^) Im J 


dv Vk^ 

-00 V — Zq 


'-cx, {v - ZqY{v - zY) \2-i j 

+ (!)' -a;(^„) + -'‘"7/_ 


(89) 


dv T]^ 


(«) 


-00 (u - ZoY 


'ik" 


S 


dv T}^ 


(*) 


2 


1-00 {v - ZoY 


6 


dv T]^ 


(s) 


-00 (v - ZoY 


Thus for < ^ ' ^ 0,0 > the typical divergence is 7 “^. 

Combining these results with pi in (|76D , we obtain the asymptotic form of the cubic 
coefficient as 7 —> O’*" 


Pi = — 
7 


Cl( 7 ) - 7 cii( 7 ) + 0 ( 7 ^) 
where ci and di are nonsingular functions of 7 defined by 


(90) 




1-00 V — Zq 


diil) = 7 




4 AK'Yzq) 


EEEi 


f 5") 

/OO ’ 

-OO [V — ^o) 


(91) 

(92) 
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this result was briefly presented without derivation in an earlier paper. |]^ When the leading 
term is nonzero at 7 = 0 (i.e. Ci(0) 7^ 0), then pi diverges like 7“^ as a result of the 
divergence in < dy'ip, k ■ ho,o >• 

The significance of this result is best understood in light of our earlier discussion of 
scaling the amplitude A{t) to obtain an amplitude equation free of singular features at small 
growth rates. With our expression for pi, the amplitude equation for the unstable mode is 


A 


A 


A + 


[ci(0) + 0(7)] 


7^ 


\A\^ + P2{7)\A\^ + 


which we rewrite using polar variables A = pe as 

P = 7 P + :^|Re(ci(0)) + 0(7)] p* + Re(p2) p*' + ■ • ■ 

e = U1- -!j|lm(ci(0)) + 0(7)1 p’ - Im(p2) p‘‘ H-. 

7 


(93) 


(94) 

(96) 


Now we introduce a scaled amplitude p{t) = following and anticipate that 5 = 1 

is required to remove the singular effect of the linear term: 


rJv ^ 

3- = r + ——[Re(ci( 0 )) + 0(7)] + 7^^“^Re(p2) - ( 96 ) 

dr 7^ 

6 = u) --[lm(ci( 0 )) + 0(7)] - 7‘^^Im(p2) -• ( 97 ) 

7 

The exponent fd should be chosen so that the nonlinear coefficients have finite limits as 
7 —>■ O’*'. The most stringent requirements on (3 arise from the coefficients in the dr/dr 
equation, and the cubic term requires 2/3 > 5 to absorb the singularity. Setting /3 = 5/2 
suffices to cancel the singularities in pi and allow the cubic nonlinearity to formally balance 
with the linear term at small 7 : 


— = r + [Re(ci( 0 )) + 0(7)] + 7®Re(p2) -\ - ( 98 ) 

dr 

0 = o; — 7 [lm(ci( 0 )) + 0(7)] — 7^°Im(p2) + ■ ■ ■. ( 99 ) 


This choice for fd assumes that Re(ci( 0 )) 7^ 0 and that the singularities at fifth (and 
higher) order are also removed by the same scaling that absorbs the cubic singularity. Typ¬ 
ically the condition Re(ci( 0 )) 7^ 0 is satisfied, but there are interesting exceptions and they 
are discussed next. The adequacy of the /3 = 5/2 scaling in controlling higher order singu¬ 
larities is studied in section m 


D. Special Cases: ci(0) = 0 

Evaluating ci( 0 ) as 7 — 0 + yields 

vr k 


ci(0) = 


4AUt7 + s0) 




(«)/ 


( 100 ) 
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where r’p(O) denotes the linear phase velocity at zero growth rate. Since the hrst factor 
(vr/c/dA^) cannot vanish, the special case Ci(0) = 0 will only arise if the species sum (which 
excludes the electrons) vanishes. By carefully adjusting the parameters of a given model, 
one can presumably arrange for exact cancellations between different terms in this sum, but 
there is little motivation to study such artihcial situations. 

However, there are three natural circumstances in which ci(0) = 0 because every term 
in the sum vanishes independently, and the cubic coefficient is clearly less singular: 

a. Inhnitely massive ions: = 0 for all s ^ e. 

b. Zero slope for the resonant ions: — 0 s ^ e. 

c. An electron-positron plasma: = 1 for positrons. 


In all three cases, the cubic coefficient has the asymptotic form 


Pi = Mi(0) + ^>(7)] 

ryO 


( 101 ) 


with di(0) given by the full formula for the second case and (ii(0) = 1/4 in 
third cases. 


In the rescaled equation (]96|), this weakened singularity requires only that 
setting {5 = 2 yields 


the hrst and 
2/3 > 4, and 


— = r + [Re(di(0)) -h 0 ( 7 )] + ■j'^Re{p 2 ) r® H- (102) 

dr 

9 = uj — [lm((ii(0)) -I- 0 ( 7 )] — 7 ®Im(p 2 ) + ■ ■ ■ ■ (103) 


In (I), we have shown that the exponent {3 = 2 suffices to absorb the singularities in ( p.02|) - 
( p.03|) to all orders for a plasma with hxed ions; this exponent corresponds to the relatively 
familiar trapping scaling for the electric held of an unstable mode. In Section |V], we provide 
a new proof of this result and also establish the same conclusion for the case (b) where each 
ion distribution is hat at the phase velocity: 'yik\vp) = 0 for all s ^ e. This will occur, for 
example, if the ions are sufficiently cold and do not populate the resonant region in velocity 
space. On the other hand, when there is an exact rehection symmetry F^iy) = Fq{—v), 
we can have an unstable mode with vp = 0 and 7 fc*^( 0 ) = 0 holds even though the ion 
distribution function is non-zero at n = 0. This latter situation arises for instabilities 
of rehection-symmetric equilibria due to real eigenvalues, e.g. a rehection-symmetric two- 
stream instability. 

Our understanding of the third example, an electron-positron plasma, is less complete. 
In Section [VB| we obtain by explicit calculation Re(p 2 ) ~ 7~® so the rescaled hfth order 
term in (lU^) remains singular unless (3 > 9/4. However, the singularity structure of the 
full expansion is more complicated and we have not determined the ehects of singularities 
at higher order. 
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IV. SINGULARITY STRUCTURE OF THE EXPANSION 


Our goal is a systematic analysis of the singularities of pj to all orders in the amplitude 
equation (pll). The detailed calculation of pj rapidly becomes prohibitively laborious, but 
the recursion relations (|57|) - (|^ determining the higher order coefficients in terms of lower 
order quantities (cf. Table I) can be analyzed to determine the properties of pj. This study 
requires an accurate estimate of the pinching singularities in the integrals found in pj. For 
this purpose we introduce an “index” which allows the divergence of a given integral to be 
assessed by a simple counting procedure. As in our discussion of the cubic coefficient, we 
assume the analytic singularities of Fq do not affect the divergence behavior; the precise 
assumption on Fq is given in (|108|) below. 


A. Definition of the index 


At every order in the expansion, the pinching singularities have a common structure. 
For n > 0, dehne 


Dn{a,v) 


1 

(n - ai)(n - 02) ■ • • (n - «„) 


(104) 


where a = (ai,..., «„) and also dehne DQ{a,v) = 1. Evaluating pj for j > 1 involves 
integrands of the form 


giv) = D^{P,vrD^{a,v)Y.{K^^^) 

S 


d<l 

dv‘> 


(105) 


with m + n > 1, m' > 1 and q > 0. The poles in g{v) may be written as 


aj = zq + i'ji'j/k j = 1 ,..., n (106) 

(3* = zl- iyCj/ k j = 1 ,..., m; (107) 

since they always lie along the vertical line Re(n) = ^^( 7 ). The specihc pole locations are 
given by the numbers Vj > 0 and Q > 0 which vary depending on the integral. In all cases, 
these numbers are independent of Fq] in particular Pj and Q are independent of 7 . 

The complex-analytic singularities of Fq are assumed to be unimportant for the amplitude 
expansions in the following sense. For any n > 1 and g > 0, we assume that 


lim 

7^0+ 



dv Dn{a, v) 


d<i 

dv^ 


< 00 


(108) 


holds for each species and arbitrary choices of Pj in ( p.06|) . For example, it is sufficient for 
each distribution function Fq^\v) to be analytic on a neighborhood of np(0), the linear phase 
velocity at zero growth rate; in this case the limit in (|108| ) can be evaluated by the Plemej 
formulas for Cauchy integrals. |]T4 


The index of g{v) 


m 


is dehned by 

Ind [Q] = m + n + q — 1, 


(109) 
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and Ind [Q] > 0 since we assnme m + n > 1. For example, for in (^), the integrand of 
r 2 ,o in (^U) and integrand of < > in (^ 0 |) , we obtain the indices 


Ind 




0 Ind 


(n - zoY 


2 


Ind 


{v - ZoY{v - Z^) 


2 . ( 110 ) 


If mn = 0, onr assnmption (|108|) ensures that / dvQ{v) is non-singular, but when mn Y 0) 
then as 7 — 0 +, the integral of Q develops a pinching singularity at the phase velocity np( 0 ). 
The index of Q estimates the worst possible divergence of J dvQ due to this singularity. 


Lemma IV.1 For G{v) in (|105D with mn Y 0, the integral of Q satisfies 


lim 7 “^ 

7—>0+ 


dvQ{v] 


< oo 


with J = Ind [Q]. However, if 

+ (k'*’)”'-'! TTp(o)) = 0, 


( 111 ) 


( 112 ) 


then the integral of Q is less singular and ( |lllj) holds with J = Ind [^] — 1. In (|112|) the 
prime on the summation indicates that the sum omits the electron species. 


Proof. 


By integrating by parts q times we can reduce to the case with g = 0. With 
g = 0 , consider the simplest possiblity m -|- n = 2 , and make a partial fraction 
expansion of the integrand 



{v - ai){v - fit) 


(2 + i>i + ^i) 



dv 


Es 4 "^ 

{v - ai) 



dv 


Es 

{v - fit) 



(113) 


The integrals on the right are non-singular with limits at 7 = 0 given by the 
Plemej formula 


lim 

7 ^ 0 + 



dv 


Es 

(n - oi) 



{v - fit) 

27rz^(«:W)-' vl^\vY0)). 

S 


(114) 


From the limit 7 — 

E. i)t7p(0)) 

(e) 


be rewritten as 7 ^ 

(EH): 


O’*" of the dispersion relation Afc( 2 ;o) = 0 , we obtain 
0 for the imaginary part. Since = —1, this can 
Es hi^^('^p( 0 )), and used to eliminate from 
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lim 

7^0+ 


- ^ Es r , v\ 


(*) 


dv 


{v - ai) 


dv 


{v - (dt) 




(115) 


This proves the lemma for m + n = 2. It is important that the limit in (115) 
is independent of the parameters z/i and Ci that locate the poles Oi and j3i for 
7 > 0. 


For integrals with m + n > 2, by expanding the integrand in partial fractions, 
they can be re-expressed in terms of integrals with m -|- n — 1 mnltiplied by a 
factor of 7 “^. A simple indnction argnment then establishes (| 111 D for general 
mn 7 ^ 0. (Recall that if mn = 0 then the integrals are nonsingnlar.) 

The general significance of the condition in (|112| ) becomes clear from the form of 
the partial fraction expansion. For example, consider the expansion for (m, n) = 
( 2 . 1 ): 


dv 




-ik/'f) 




-oo {v - ai){v - a 2 ){v - I3i) {2 + vi + (i) J-oo {v - ai){v - 02 ) 


(116) 


-ik/'^y 


[2 + i>i + Cl) (2 + V 2 + Cl) 


dv 


{V - Q!2) 


dv 


{v - /?r) 


The second term is the dominant singnlarity 7 “^ and at 7 = 0 the bracketed 
integrals yield the same factor obtained previonsly in (p.l5|) . When this factor is 
zero, the singnlarity of the {m,n) = ( 2 , 1 ) integral is rednced from 2 to 1 . 

In the partial fraction expansion of the general case, the most singnlar terms 
that diverge according to (| 111 |) will always be proportional to the difference of 
two non-singnlar integrals: 



dv 


Es # 

(n - aj) 



dv 


Es # 

(v - l3p) 


(117) 


whose limit at 7 = 0 is given by ( p.l5| ). Thns when this limit is zero, the 
singnlarity of the term mnst drop by 1. This proves the second part of the 
lemma. 

□ 


The exceptional sitnation defined by the condition in (|112|) is a generalization of the 
featnre noted previonsly in onr analysis of the cnbic coefficient. Setting m' = 3 in (|112|) 
yields ci(0) = 0 from ( p.0CI|) ; the special circnmstance that corresponded to a less singnlar 
integral at third order. 

Onr application of Lemma [lV.l| to the recnrsion relations for reqnires generalization 
of the index to allow for snms of fnnctions with well-defined indices and prodncts of Q 
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with singular functions of 7 . In each case, the generalized index is dehned so that Eq. (|111|) 
remains true, i.e. the index for the composite function gives the maxima/possible divergence 
of its integral. 

First, if Qiiy) and ^ 2 ( 1 ^) have indices satisfying Ind [^ 1 ] > Ind [^ 2 ] then we define the 
index of the sum to be the larger index: 

Ind [ 6/1 + 6 / 2 ] = Ind [e/i], (118) 

Clearly, Ind [Qi] gives the maximal possible divergence of / dv {Qi + ^2)- Secondly, if q{pf) 
is a function of 7 with the asymptotic behavior ^( 7 ) ~ 7 “^ as 7 —O’*", then we dehne the 
index of ^( 7 ) Qiy) to be 


Ind [q Q] = Ind [G] + m 


(119) 


Estimates of the form (|111|) still hold for qQ with J = Ind [q G\. This completes the dehnition 
of the index. 

By applying (|118| ) and ( |119|) , the indices of and niay be determined. For 
example, from and we have the indices 


Ind 


and from (|77D and (|79|) 


Ind 






= 1 


Ind 


= 2 


Ind 






= 1, 


( 120 ) 


= 2 . 


( 121 ) 


The index for h^ 2 fl is correct provided r 2 ,o does not diverge more strongly than 7 “^. We 
can conhrm this immediately from the index in (| 110 |) ; in fact, r 2 ,o is non-singular, as noted 
already in Section [III C| . In section jlV Bj , these indices are applied to determine the singu¬ 
larity of pi. We stress again that when applied to a composite function G{v) the estimate in 


( 111 ) does not necessarily determine the true singularity, but only an upper bound on the 
possible divergence of the integral. 

There are several immediate consequences of the index dehnition worth stating. First, 
complex conjugation doesn’t alter the index, Ind [G] = Ind [G*]', secondly, if G{v) has a 
well-dehned index, then dividing G by (n — a) or {v — (3*) simply increases the index of G 
by one: 


Ind [G/{v - a)] = Ind [G/{v - jS*)] = Ind [G] + 1. 


( 122 ) 


Here a,/? are dehned as in ( [106] ) - ( jlOTj) and thereby increase the index of every term in G 
by one; this implies (| 122 [ ). Finally, from (|105|) , diherentiating G also raises the index by one, 
Ind [d^G] = 1 + Ind [G]. 


B. Analysis of pi using the index 

The above relations along with recursion relations determine the indices of Im,j and 
to all orders. This information permits an estimate of the singularities of the coefficients 
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Pj for j > 1. In practice, this must be done recursively following the organization of the 
amplitude expansions in Table I. The procedure is illustrated here for the hrst level of Table 
I; more general results on the index to all orders are provided in Section |1V C| and Section 

0 

Table I begins with -0, ho,o and ^ 2,0 whose indices were evaluated in (|11(]| ) and (|121D . 
From this information, we can estimate the integrals required to evaluate pi in (^). The 
integrand of < dy'ip, k ■ 'ip* > has index 2, and the integrands of < dv'ip, n ■ ho,o > and 
< K ■ h 2 fi > have index 4. In each case, the factor increases the index by two 

from the values obtained in (IIU) and (1^). Now Lemma |IV.lj tells us that < dy'ip, n-pj* > 
cannot diverge more strongly than 7 “^, and < k ■ > and < k ■ h 2 fl > cannot 

diverge more strongly than 7 ““^. 

These conclusions imply that the singularity of pi cannot be worse than 7 “^^; this was 
found in Section [III Q by an explicit calculation of pi. The explicit evaluation, of course, 
provided the important additional result that the 7 “"^ singularity actually occurs unless one 
of the exceptional cases corresponding to Ci(0) = 0 is considered. The foregoing index 
analysis is easier but yields less complete information. 


C. The Main Result 


The main result on the singularity structure of the amplitude expansions can now be 
proved. 

Theorem IV. 1 For j > 1, the coefficients in the expansion of the amplitude equation (|54D 
satisfy 


lim 7 ®-^ ^ \pj\ < 00 . (123) 

7 ^ 0 + 

Let Jm,j = (2m + 5j — 3) + 45^,0 + 56m, 1 , then for j > 0, m > 0 and m' > 0, the indices of 
I'r and h'-r obey 


Ind 

S 

Ind 

_ s 

and the integrals in (^) and ([75|) satisfy 


— Jm,j- 


— Jm,j T 1, 


+ l Ip .1 ^ 
I-*- m,j I ^ 


lim ^m,l 

7—>0+ 

(«)/ 


lim 7 

7^0+ 

7 dyf ^ K * hm,j ^ 


< 00 . 


In addition, integrals over If I (v) satisfy 


(124) 

(125) 


(126) 

(127) 


lim -/‘o+n 1 

7^0+ 




< 00 


(128) 


for n > 1 where Dn{a, v) is defined in ( |104|) with upper half-plane poles a = (cti, 02 ,..., af) 
of the form (|106|) but otherwise arbitrary. 
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The relations in (|124|) - (|128|) are of interest chiefly for their utility in proving ( p.23|) . We 
are able to prove (|124|) -( |I^) as upper bounds, but we expect them to hold as equalities in 
most cases. For m = 1, the bounds in (|126D and ( |127|) give stronger estimates on the integrals 
than one would obtain simply from the index of the integrand calculated using (|124|) - (|125|) . 
Similarly, the bounds in (|128|) are needed because they provide sharper estimates for the 


m = 1 case than are available using the index information (|124|) - (|125|) alone. 

Even though hij and have the same index, it is very important to recognize that 
the conclusion in (|127D may not hold for < dv'ip,K ■ h\j >. The most singular terms in 
< dv'ipyH ■ hij > are not in general the most singular terms in < >; for < 

■ hlj >, if we wish to avoid explicitly evaluating the integral, our only means of 
estimating the divergence is by calculating the index of the integrand from 
Proof . 


1. The proof is by induction using the organization of Table I. In Section IV A 


IV B|, the relations i 


m 


(|127|) have been explicitly verihed for Iq o, I- 


2 , 0 ; 


^ 0 , 0 ) ^ 2 , 0 ) To^o, r 2 ,o, and pi. This proves the theorem for the first level of 
Table I, but does not establish the validity of (|128|) for the initial case j = 0 
which enters in the second level of Table I. By substituting the recursion 
relation (^) for Ji^o into ( 128 ) we obtain 


7 


3+n 


dv DJa,v)Y,(f^''Y li%’) = 


(129) 


7 


3+n 


kj J-c 


dvDn{a,v)d^ [k- [ho,o - ^2,0 + -^'^*'^2,0] 


+7 < dviJ, K • [ho,o - ^ 2,0 + 2 ^*r 2 ,o] > J dv v)ip{v). 

On the right hand side, the index of the first integrand is {n + 3) so the 
factor of 7 "+^ ensures a hnite limit. In the second term, the integral over 
Dn'il’ is non-singular since all poles are in the upper half-plane, but the 
coefficient < dy'ild, n ■ [ho,o ~ ^ 2,0 + \'^*'^ 2 fi] > diverges like 7 “*^ (cf. Section 
111 C). Thus for n > 1 the second term will also have a finite limit; this 


verifies (p. 28 |) for Ji^o- 


2 . 


The relations 

and pj using the recursion relations (0) 


( p.28|) are extended by induction to all coefficients 1 ^,^, 

and m - o- Assume 
that (|123|) - (|128|) are valid down to some arbitrary level of Table I, and 
consider what the recursion relations imply for the coefficients evaluated at 
the next level. We use to denote subscripts of coefficients such as 

and pj^ that are to be evaluated from lower order quantities 
assumed to satisfy (|123|) - ( p.28|) . 

The induction argument has three parts. First, we show that the properties 
( p.25|) - ( [I27|) must hold for hmi, 
namely the index identity (|124|) 


- ( [I27|) must hold for hmiji and if Imiji satisfies two relations. 


Ind 


S 








(130) 


26 













































and the irti = 1 estimate (| 12 ^) 


lim 7 

7—*0+ 


Jl,j, +n—1 


dvDn{a, V 




j(«) / 

-'ijd 


< 00 . 


(131) 


In the second part of the proof, we verify that (|130|) - (|131|) do hold for 


-‘mi Jl- 

imply 


Finally, in the third part of the proof, we verify that 


(E) 


for pj^. It is important to note from Table I that the quantities 
miji depend on pj for j < ji, but are independent of pj^ so the reasoning 
is not circular. 


h 


(Part 1): Suppose that (|130|) - (|131|) hold, then the three relations (|125| ) 
- ( p.27| ) are easily obtained as follows. First, consider ( |126D which holds 
trivially for mi = 0 since Fqji = 0. For mi > 2, (|126|) follows from (|130| ) 
and the identity (^). For mi = 1, we have to allow for the fact that 
Afc(ziji) ~ 7 as 7 —>• O’*", and apply (|131|) to estimate the integral in (|7^). 
Next ( |125| ) follows by applying ( |13U| ) and (|12ti| ) to the identity (]75|). Each 
term in (|73[) separately has a maximum index of {Jmiji + 1) so this must 
bound the index of the sum. 


Finally consider ( [127|) . For mi 7 ^ 1, this relation follows from (|125|) , and 
for mi = 1 we evaluate (|127|) using the identity in ([T3|) : 


K ■ hij7 > = < d^ip, K • hjJiv - ZijJ > (132) 

< dyTp, ■ Pk/(v - ZijJ > 


The hrst term on the right is non-singular by virtue of ( 131 ), and the second 
term is also non-singular by virtue of 


4. (Part 2); Thus the crux of the matter is to verify (|130|) and (|131|) from the 


recursion relations for i-®- 


from (| 66 D - (^). In this step, we assume 
( p.23|) - ( p.28|) apply to the quantities appearing on the right hand side of 
these recursion relations, then the properties in ( |13CI| ) - ( p.31|) are verihed 
term by term. 

Verification of im - the index of 

The index of each term on the right hand side of the recursion relations 
for Imljii'o) can be evaluated by applying ( |123|) - (|127| ); this exercise shows 
that all of these terms have an index less than or equal to Jmi,ji which 
establishes (p.30| ) for A few examples from (| 68 D, the recursion relation 

for (mi,ji) = ( 1 , ji), illustrate these index calculations; note that in these 
examples Jij^ = 5ji -|- 4. 


The terms in 

M. 


that depend on pj have the form [(2 -|- n)pj^-n + (1 + 


from (|125|) we have 


Ind 




^ ^l,n “t“ 5, 


(133) 
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and from (|123|) the singularity of Pj^-n is determined, hence 


Ind 




< 5(ji-n) - 1 +(5n + 5) = 5ji+4 (134) 


which is consistent with ( |130|) for = (1, ji). 

(b) Next consider the term ■ dy in (^); from the dehnition (^) 
we have that the index of V_lk ■ dy is equal to the largest index 
obtained from the two terms n ■ dyhoj-^ and < dy'ipjK ■ >. 

These terms have individual indices 


Ind 




Ind 




< 5ji + 3 

< 5ji + 4 


(135) 

(136) 


where the first index 
from (|127|) . Hence 


comes from 


and the second follows 


Ind 




< 6ii + 4 


(137) 


which is consistent with ( |130|) for (mi,ji) = (1, ji). 

(c) A more subtle example is the term V± k- dyh\ ^T2ji-n-i in 
expands to 


which 


'PK • dy ^ ■ dy h^ j^T'2Ji—n—l 

+r2ji-n-l V’(n) < dy'ljj, K ■ > . 

Now the index of each term can be evaluated to give 


(138) 


Ind 


Ind 




thus 


Ind 




jl-n-l 


< 5i, + 2 

< 5ii + 4; 

< 5ii + 4 


(139) 

(140) 

(141) 


which is consistent with (1^) for (mi,ji) = (1, ji). In this last ex¬ 
ample, ( |139|) is obtained from ( |125|) and ( [I26|) ; the second result 


differs from (|136|) since (|127|) cannot be applied and the singularity of 
< ^y^jJ, K ■ h\^> must be estimated from the index in ( |125| ). 

In a similar way the index of each term appearing on the right in the 
recursion relations (|66D - (fflD for /^^^(n) can be estimated with the hnal 
result in (|130|) . 
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5. Verification o/O); 

We integrate the recursion relation for and examine the resulting 

expression to establish ( |131|) . When (| 68 [), the recursion for is inserted 


into (|131|) we obtain 


7 


Jl,ii +n—1 


dvDr,{a,v)J2{K^^Y liiiv 


ii-1 


[(2 + 0PA-i + (l + 0Ptd / dvD4a,v)J^{i.<‘Yhi%') 

1 n ^ J—OO c 


+7' 


z=o 


5jl+n+3 




(142) 


+^5ii+n+3 ( 7 ) < K ■ G > 

\k; 


dvD^{a,v)Y,{K^^Y 


where G{y) = hoji — h 2 j^ + ■ ■ ■ denotes the collection of terms enclosed 
by brackets in (|68D . The relations ( p.24|) - (|127|) hold for these terms by 
assumption, and it is straightforward to verify that Ind [G] < 5ji + 2; 
this immediately implies that the terms in (|142|) involving G are finite as 
7 —> 0^. For the remaining terms in (|142|) we have Pj^-i ~ from 


(p.23|), hence we need to prove 


lim 7 '^'+^+’^ l ~ dvDn{a,v)J2YY h["}{v) 


< OO 


(n > 1) (143) 


to establish (|131|) . The formula in (|125|) bounds the index of the integrand 
in (|143| ) by only (5/ + n + 5) which is not sufficient. Instead, we evaluate 
using the identity (|73[) 




n r , r,, 

- / dvDn{a,v) - - -r- 

kj J -00 [V — Zifi 

-Ti,, / dv D^ia, v) -A. (144) 

7-00 (r> - Zifi 


The integral in the second term is non-singular and Ti^^ ~ 7 ® from 


thus the desired estimate in ( p.43|) holds for the second term. It remains to 
verify 


lim 

7 ^ 0 + 


dv Dn(a, V 




- Y 


< 00 


(n > 1) (145) 


for the first term in (|144D ; this follows from 

Dn/{v - Zifi. 

The proof of ( p.43| ) establishes ( |131|) for G ^ 7 . 


with j I and 
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6 . (Part 3): Verification of (|123|) for 

It is straightforward to check that application of (|125| ) - ( p.27|) to the right 
hand side of the recursion relations (|F7| ) - (|59|) yields the estimate in ( |123| ). 
A few examples from (|5^) - (|5^ suffice to illustrate how this is carried out. 

(a) Consider the term < K-hoji-i > in (^); from (p.27|) the singularity 
of this integral is at most 


< K ■ ho,,1-1 > 




t). 


(146) 


which is consistent with ( |123|) . 

(b) Other terms will involve products of integrals; for example, r 2 ,ji-i < 
> in (^. Using (| 8 ^) and (|126|) we determine that the 
singularity of this product is at most 








(147) 


which again is consistent with (|123|) . 

(c) A hnal example is the product V 2 ,i < n ■ > with 0 < / < 

and 


ji — 2 from ( 1 ^ ; applying 
singularity of 


^2,1 < dyi!,K- > 


to this shows a maximum 


^1 \ l + /-I \ _i_2 

t it 


5ji-l 




(148) 


which is consistent with 
This completes the proof. 


□ 


D. Implications of the main result 


Theorem m allows us to evaluate the 7 ^ O’*" limit of the theory. Rewriting the 
amplitude equation (|^) using polar variables A(t) = p(t) exp{—i9(t)), and the scaling 


from Section mg, yields 


= r 


dr 
dr 

de 

-dt=^-^ 


00 

El 


7 




CXJ 

El 


7 


®/^r(7f) 

(149) 

“^Re(pj) 

(150) 

fi^~Am.{pj) r‘^fi 

(161) 


where r = yf. The result ( [123|) on the divergence of the coefficients pj ensures that the 
terms in (p.5(]| ) - (|151|) are now hnite as 7 —> O’*". Introducing the asymptotic coefficients 
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(162) 


Cj{0) = lim 

7—>■0+ 


7^^ Vi 


(j > 1 ) 


the amplitude equation for a weakly unstable wave becomes 

dr 


dr 


= r + [R-6(7(0)) + 0 ( 7 )] 


(153) 


i=i 


Unless Re(cj(0)) = 0 in each term, the scaling in (|149|) is required to obtain a hnite theory. 
The circumstances that could force ci(0) = 0 were discussed in Section P^II C| , and will be 
analyzed more completely in Section 

The electric field is obtained from Poisson’s equation (imk) Emk(t) = J2s 
using the form of in (|4^) 

iA{l + aTi{a))/k m = 1 

(154) 

iA^ Tm{cr)/mk m> 1. 

The small growth rate behavior of rm(cr) can be inferred from ( |126| ) and the expansion 
rm(o') = Using we hnd an expansion of T^ with a non-singular 

limit 


Emkit) 


2m-2+5Sm,ir — r 

I ^ m / , 1 ^ m,j ' 


27 


j=0 


and this motivates our definition of TIL 


T^ = hm 72 — 


7 ^ 0 + 


j=0 


where 


bm,j{0) = lim 

7^0+ *- 


I ^ ' 




(155) 


(166) 


(167) 


With this notation, the asymptotic form of the Fourier expansion E{x,t) = 
J2mEmkit) exp{imkx) of the held can be evaluated 


E{x, t) = 


,^ 5/2 


r(r) (1 + r;r(r)‘‘ + 0(7)) 

(r;;.r(r)’" + 0 ( 7 )) 


(168) 


m=2 


m 


+ cc. 


The asymptotic form may be rewritten more simply as 


iE{x, t) 


7 


5/2 


r{r) 

k 


(1 + r7(r)^) -f cc 


+ o{^y, 


(159) 


this describes an electric held comprised of a single Fourier component whose coefficient 
scales like 7 ®/^ and evolves nonlinearly according to the asymptotic amplitude equation 
The 7 ^/^ scaling yields a dramatically smaller held than would be expected from the 


trapping scaling 7 ^ characteristic of a model with inhnitely massive ions. 

An analogous discussion of the asymptotic distribution function f{x,v, t) is substantially 
more complicated because of singularities in v that emerge as 7 —> O’*". This issue is briehy 
treated in Section 0 , but we defer a full analysis to a companion paper. 
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V. SPECIAL CASES 


Section ETTg identified three settings for which Ci(0) = 0 and the 7 ^ singularity of pi is 
replaced by a 7 “^ divergence: 

a. Inhnitely massive ions: = 0 for all s 7 ^ e, 

b. Zero slope for the resonant ions: Vk\vp{0)) = 0 for all s 7 ^ e, 

c. An electron-positron plasma: = 1 for positrons {s = p). 

From the partial fraction expansion in Lemma we found the general condition required 
to reduce the singularity of a standard integral, 


/ oo roo _ 77^*^ 

dvg{v)= dv Dm{g,v)* Dn{a,v) Y, r, [ , (mn > 0 ) ( 160 ) 

-OO J—OC c UV^ 


from 7 " '?+! to 7 namely, 


-)bK(o))=0 


(161) 


where the primed sum omits the electrons. The specihc instance of ( |161|) that yields ci(0) = 0 
corresponds to m! = 3. 

The validity of this general condition varies between the three settings listed above. For 
both (a) and (b), (|161|) holds quite generally for all integers m' > 0 ; however for the third 
setting (c), the condition reduces to 

|(-l)'”' + l],7W0)) = 0. (162) 


and this will be generally satished only for m' odd. 


A. Cases (a) and (b) 


It is relatively simple to adapt Theorem |1V.1| to accomodate the hrst two settings where 


( p.61|) holds uniformly. The following result restates Theorem IV. 1 of (I) under more general 
assumptions (note that the dehnition of the index in this paper differs from (I)). [|l| 


Theorem V.l Assume that (|161| ) holds for all m' > 0; then for j > 1, the coefficients in 
the expansion of the amplitude equation (|5l) satisfy 


lim ^ |pj| < cx). 
7 ^ 0 + 


(163) 


Let Jm,j = (2m - 1 - 5j — 3) -|- 45^,0 + 56m, 1 , then for j > 0, m > 0 and m' > 0, the indices of 
Im]j and hm]j obey 
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Ind 

_ S 

^ Jm,j j 

(164) 

Ind 

s 

^ Jm,j 1 j ~ 


(165) 

(^) and ([75|) satisfy 




lim ^ \Tm,j 

7 ^ 0 + 

< 00 

(166) 

lim 

- . r\4- 


A 

S 

V 

< 00 . 

(167) 


The proof is simpler since 
for special estimates such as 

Proof. 


(|167|) follow directly from (|164|) and and there is no need 
in the m = 1 case. 


1. The induction proof follows that for Theorem |1V.1| , and we only highlight 
the new features brought in by (|161|) . The recursion relations for Im,j are 
hnite sums; each term of which has the form ^( 7 ) Q{v) where Ind [^(n)] > 1 
and the coefficient ^( 7 ) may be singular as 7 —O’*". According to Lemma 
IV. 1| , with (|161|) in force, if J = Ind then 7'^“^g(7) / dvQ{v) 

remains finite as 7 —> O’*". 


2. We have already shown that (|161|) implies (|163|) for j = 1 so the theorem 
has been proved for the first level of Table I. As before, we assume that 
( p.63|) - (|167|) are valid down to some arbitrary level of Table I, and consider 
what the recursion relations imply for the coefficients evaluated at the next 
level. Again, we use to denote subscripts of coefficients such as 

Imiji, ^miji and pj^ that are to be evaluated from lower order quantities 
assumed to satisfy (|163|) - ( |167|) . Following the organization of the previous 


proof, we first check that (|164|) for I^lji implies (|165|) - (|167|) for 


Tmiji, and < dyijj, k- >, respectively. Then (|164|) is verified from the 

recursion relation for and finally (|163|) is shown to follow from (|165| ) 


- (m)- 

3 . (Part 1): Assume that ( p.64|) holds for /, 


(d . 

"U j'l' 


Ind 


S 


{k 




r('*) 


< J, 




3l 


-'mi,!) 


( 168 ) 


then the properties (|165|) - (|1671) are easily obtained as follows. First, the 
integrand of Fmiji in (|72[) has index 


Ind 



< J, 




Jl 




+ 1 


(169) 


from (|164|) so the expression 
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ryJm-^,31 


—jl—<5rr 


dv 




V — Z; 


(170) 




remains finite as 7 —O’*"; this proves (|166|) since the < 5 ^ 1,1 in the exponent 
compensates for the added singularity of in (^). Next, (|165|) fol¬ 

lows by applying ( |168D and (|166|) to the identity dZD- Finally, the integrand 
in (|167|) has index 


Ind 




is) 


K 


is)uis) 


— Jmiji jl 




+ 3 


(171) 


4, 


from ( |165|) so the bound in (|167D follows from ( p.61|) and Lemma |1V.1] . 
(Part 2): Verification of (|168|) - the index of 


The consequences of (1^) for the index calculation are incorporated in 
( p.63| ) - ( [166| ). The index of each term on the right hand side of the recur¬ 
sion relations for Imljfiv) can be evaluated by applying these conclusions 
regarding the lower order coefficients; this is done exactly as in the proof of 
( p.3CI|) in Theorem [1V.1| . The bound in ( p. 68 |) follows as before. 

(Part 3): Verification of (|163|) for pj^ 

It is straightforward to check that application of (|165| ) - ( p.67| ) to the right 
hand side of the recursion relations (^) - (|5^) yields the estimate in 
just as in the proof of Theorem [IV. 1|. 


□ 


In light of the discussion in Section jlV Dj , the theorem tells us that the scaling in (|149|) 
is replaced by 


p(t)= 7 V( 7 t) (172) 

for the special cases (a) and (b). This modification in turn implies an electric field charac¬ 
terized by the “trapping scaling” Ek ~ 7 ^; the discussion in (I) provides detailed expressions 
for the asymptotic form of E{x,t) and f^'^\x,v,t) in the limit 7 —> O’*" when the ions are 
fixed. 0 


B. An electron-positron plasma 


The result pj ~ for the two special cases just considered is not correct for the 

electron-positron system with nd’i = —= 1 . Although pi ~ 7 “^ does hold, we find at 
fifth order p 2 ~ 7 “® instead of p 2 ~ 7 “^. Thus the trapping scaling of the amplitude in 
( p.72|) is not sufficient to obtain a finite expansion. We summarize the evaluation of p 2 from 
(13) " (HD 5 tint leave the comprehensive analysis of this case as an open problem. 

The exact recursion relation for p 2 is 
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( 173 ) 


-ikp2 = - < dy^j, K ■ (/io,i - /i2,i) > “2 ^ > ^ 2,1 

- < dy^jj, K ■ hofl > ri^o+ < dy^jj, K ■ /i2,o > r^o 
- ji < dy-^, K ■ hl o > 12,0 + < dyi), K ■ h^Q > 


2,0 


< dy-lp, K ■ /l3,o > \ , 


and, in the bracketed terms, r 2 ,o and Fso are non-singular so Theorem |1V.1| ensures the 
divergence of these terms cannot exceed 7 “'^; hence they are subdominant: 


-ikp2 = - < dy-ip, K ■ (ho,i - h2,i) > “2 ^ > ^ 2,1 

— < dy'ljj, K ■ /lo,o > ri,o-|- < dy'ljj, K ■ ^2,0 > Q -|- 0(7 ^). 


(174) 


Here consider the two terms involving Ti q. In the general case. Theorem m gives ri,o ~ 
7 “®, but this drops to 7 “^ when we take (|161|) into account for m' odd. The integral 
< dy'il}, K ■ ^ 2,0 > is non-singular so < dy'ljj, n ■ h 2 ,o > ~ 7 “^^ and we neglect this term 

also. Next < dy'ljj, k ■ ho,o > is written out to hnd 


< dy'ljj, K ■ ho,0 >= ^ < dy'lp, ■ dy{'lp* - Jp) >', 


(175) 


this integrand has index equal to 3, but contains an odd power of n so (|161|) holds and the 
divergence cannot exceed 7 “^. Overall, < dy'lp, k ■ ho,o > ri,o has a maximum singularity of 
7 “”^. A careful examination of the hrst term in (|174|) yields the same estimate < dy'lp, k ■ 
(ho,i - h2,i) >~ 7"^- 

The dominant singularity in p 2 arises from the remaining term 


P2 = 


^r2,i 


< dy'lp, K-'lp* > + 0(7 ^), 


(176) 


and from (|85|) and (|115|) , we obtain 


< dy'lp, K ■ 'Ip* > = 


'k' 


AP{vp + iO) 


The evaluation of r 2 ,i starts from (|7^ 

—i/2k 


r2,i = 


E 


r~ 4a'(!') 

dv 


+ 0 ( 7 ) 


(177) 


A2fc(2:2,l) s V — Z2^l^ 

and the recursion relation (^) for / 2 ,i 

i d 


(178) 


X U ( 1 1 ^ 

h,i{v) = —2pi h2,o + ^ ■ |hi,o + V’ri,o — ^3,0 + -'0*r3,o -|- -ho,or2,o| • 


(179) 


Note that A 2 fc(^ 2 ,i) = 3/4 -|- Oj'j) from (^) so any singularities arise from the integral 
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■oo 


(180) 


dv -^-= —Zpi 


^A‘Uv) 


V - 2 : 2,1 


/ dv ■ 

1-00 V — 2 : 2,1 

dv 


, 12 E''"’ A« + i’‘-A.o-hii + A‘Afi + WoAo 

fv J —00 yu — ^21) s L o z 


The strongest singularity here occurs in the integral over and all other terms in (|180|) 
are subdominant; thus we can focus attention on the expression 


dv 


— ' dv 


V - 2 ) 2,1 


k J-c 


{v - Z 2 ,iy 


+ 0(7 


- 4 \ 


(181) 


When is evaluated using ([73| ) we obtain 

, Esl2,i{v) 1 , 

/ dv -^-= -^ / dv 




V - Z 2 ,l J-00 (v - Z 2 ,l)‘^{v - Zifi) 


+ 0(7 


- 4 \ 


(182) 


Now the recursion relation (^) for is applied, and examination of the various terms 
shows that the dominant singularity arises from the occurrence of ip* in /io,o so that ( 178 ) 
can be reduced to 


r2,i = 




±^) r 

37P ) J-00 [V - Z2,iY[v - Zifl) 


(183) 


Finally a partial fraction expansion of the remaining integral yields a 7 ^ singularity 


-00 ^ {v - Z2,iY{v - 2:1,0) 7® 


fc^4^^(np(0)) + O(7) 


(184) 


and leads to the asymptotic form of p 2 from (|176|) 


P2 = 


ryO 


i7i^k\ privp) 


(p)i 


48 J Ap{vp + iO) 


+ 0{^) 


+ 0(7- 


(185) 


Since ImA),(np + i0) is typically non-zero, this shows that Re(p2) ~ and implies that 
the trapping scaling (|172|) leaves a residual singularity at fifth order. With p(t) = 7^r(7t) 
the rescaled equation becomes 

^ = r + 7^^“^Re(pi) + 7^^“^Re(p2) r® + • ■ ■ ( 186 ) 

dr 

so /3 > 9/4 is required to overcome the 7“® divergence in p 2 . 


VI. DISCUSSION 

Our analysis of an unstable electrostatic wave expresses the Vlasov distribution function 
f{x,v,t) in terms of the wave amplitude, i.e. 
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f{x,v,t)= A{t)e^'^=^^/J{v) + cc 


+ 


Amkx 


H^k{,v,A,A* 


(187) 


where A{t) evolves according to an amplitude equation, 

A{t) = Ap{\A^\). 


(188) 


In the small growth rate limit, the eigenvalue of the unstable mode approaches the imaginary 
axis and, at 7 = 0 , merges with a neutrally stable continuous spectrum on the axis. iffl 
This “interaction” between the unstable mode and the neutrally stable continuum leads to 
singular features in the asymptotic form of the amplitude equation (|188|) 


r| 


These singularities are physically signihcant since they reflect the very strong nonlinear 
interaction between the wave and the resonant particles at the linear phase velocity. Re¬ 
markably, the quantitative form of the singularities can be consistently interpreted as hxing 
overall scalings of the nonlinear solution with the linear growth rate in the asymptotic limit 
7 —>■ O’*". This interpretation is achieved through a singular transformation |R(t)| = r{'yt) 
of the amplitude in which the exponent (3 is hxed by the criterion that the rescaled dynamics. 


dr Re(p(7^^r(r)^)) 

—— ; 

dr 7 


(189) 


is free of singularities as 7 —O’*". Our conclusions are summarized in Table II; with the 
notable exception of the electron-positron system, we can rigorously establish an exact value 
oi f3 = 5/2 in the typical instability and [3 = 2 for certain exceptional models. In the electron- 
positron case, the we have hxed only a range of possible (3 values from the explicit calculation 
oi p{\A\A through fourth order in |A|; the complexity of the recursion relations in this case 
prevents the iterative determination of the dominant singularities at each order. 

The shift from /? = 2to/d = 5/2as the true asymptotic scaling due to hnite mass ions 
at the phase velocity of the wave is our most signihcant conclusion regarding the physical 
character of the nonlinear wave. Since previous studies of the asymptotic scaling of single 
wave instabilities have focussed primarily on plasma waves and assumed hxed ions, this 
shift was not anticipated and was discovered purely as a consequence of the singularities 
in the amplitude equation ( |188|) . Numerical simulations of unstable ion-acoustic waves or 
unstable plasma waves in systems with hnite mass ratios would be of considerable interest in 
this connection. In particular, for plasma waves one expects a crossover in the asymptotic 
scaling from the trapping scaling (/? = 2) to /? = 5/2 at sufficiently small growth rates. 
Predictions of this crossover have been made elsewhere based on the competition between 
the 7 “^ and 7 “^ singularities in the cubic coefficient. 

In this paper we have not analyzed the asymptotic form of /(x, n, t), but it is clear from 
the expansions for H{x, v, A, A*) that this is an interesting question. For example, consider 
the leading coefficient in ho{v) from ([77|), 


h ( \ ^ ^ 

= -l^Tv 


■ rjk 


{v - zo){v - z^) 


(190) 


there are poles in the complex-velocity plane at Zq and Zq that approach the real velocity 
axis at Vp in the limit 7 —> O’*". Thus the coefficient becomes a singular function at the phase 
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velocity when 7 = 0 ; on the other hand, for v 7 ^ Vp, the coefficient is well-behaved at small 


growth rates. In the distribution function (|1871) this particular coefficient is multiplied by 
so that after rescaling the amplitude we have a term of the form Now the 

emerging singularity in velocity space is weighted by a factor of 7 ^^, and one would like to 
understand whether the resulting product is hnite or singular at 7 = 0. A detailed analysis 
of this problem will be presented in a forthcoming paper, but the primary conclusion can 
be simply summarized. [|^ Inside a neighborhood of width and centered on n = Vp, 

the singularities in the coefficients of H are balanced by the factors of 7 introduced by 
the rescaling of |A| and one can obtain hnite, non-zero limits for the weighted coefficients. 
Outside this small neighborhood of the phase velocity, the factors of 7 force the weighted 
coefficients to vanish at 7 = 0 . 

The connection between coefficient singularities and nonlinear scalings that is studied in 
this paper is relevant in other settings. As mentioned above, a key feature of our instabilities 
is the emergence of the critical eigenvalues from a neutrally stable continuous spectrum. It 
has been appreciated recently that this characteristic is shared by a diverse set of additional 
examples which include the onset of linear instability in huid shear flows m- E 


insta- 

synchronization of large systems of coupled oscillators, 
^ and bifurcations in “mean held” descriptions of the dynamics of bubble clouds in huids 


bilities of solitary waves |21 


In the phase-dynamical models of coupled oscillators, the amplitude equations of the 
weakly unstable modes have been analyzed in the 7 —> O’*" limit. In these studies, the 


coefficients can be singular or non-singular depending on how the oscillators are coupled, 
but when singularities occur in the limit of small growth they can be absorbed by singular 
rescalings of the amplitude in a procedure analogous to the method used in this paper. The 
specihc values of the scaling exponents are different, but yield predictions for the emergence 
of the nonlinear synchronized state. 

It seems likely that the analysis of singular amplitude equations provides a general means 
to determine scaling exponents for weakly unstable systems. The exponents in turn provide a 
quantitative measure of the strong nonlinear interaction between the neutrally stable modes 
and the weakly unstable modes. 
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TABLES 


TABLE 1. Order of calculation of hm,j{v) and pj from ip{v). The flow of calculation of the 
hm,j{v) is indicated by moving downward. From 'ip{v), /io,o and /i 2 ,o can be calculated and then 
Pi determined; hi^ and are calculated next from {pi, ho,o, ^ 2 , 0 } and then /io,i and /i 2 ,i can be 
evaluated. This then determines p 2 , and so forth. For ji > 2, pj^ requires prior calculation of hm,j 
for 0 < m < ji + 1 and 0<j<ji-m + l- 2(5m,o + 



m = 0 

m = 1 

m = 2 

m = 3 

m = 4 

m = 5 

m = 6 

Po 

V>(u) 

Pi 

^ 0,0 

. 

h 2,0 





P2 

^ 0,1 

^ 1,0 

h 2 ,i 

h3,0 




P3 

ho ,2 

^ 1,1 

h 2,2 

h3,i 

hA ,0 



PA 

ho ,3 

hi ,2 

h2,3 

h3,2 

h/i,i 

h5,0 


P5 

ho,A 

hi,3 

^ 2,4 

h3,3 

^ 4,2 

^ 5,1 

ho,o 




TABLE II. 
(b) and (c) are 

Summary of the exponent /3 obtained from the amplitude equation. The cases (a), 
the exceptional cases discussed in Section ^ 

System 

P 

Comment 

ci(0) /O 

5 

2 

generic result for multiple mobile species 

(a) ci(0) = 0 

2 

infinitely massive ions 

(b) ci(0) = 0 

2 

(vp) = 0 for each ion species 

(c) ci(0) = 0 


an electron-positron plasma 
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